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Abstract: An invariance principle for Azema martingales is presented as well as a new device 
to construct solutions of Emery's structure equations. 

1. Introduction 

In Meyer raised the problem of finding "normal" martingales, i.e. martingales Mt 
on a filtered probability space {0,,J^,J^t,P) satisfying {M,M)t = t, which would enjoy 
the chaotic representation property. This property, which is stronger than the previsible 
representation property, requires the direct sum of the chaoses of M to be equal to the 
whole space L'^{Q). The only known examples of such martingales were the Brownian 
motion and the compensated Poisson process until Emery found in [3] a way to construct a 
whole family of examples, by introducing the beautiful device of structure equations. Among 
these examples, some of them, which are self-similar, appeared to be the generalization of 
two martingales considered earlier, for different purposes. One had been considered by 
Azema in pP, and the other by Parthasarathy in an unpublished paper ilSj. The whole 
class of these self-similar martingales is now usually called Azema martingales. However, 
their behavior remains a little mysterious, especially, in the neigbourhood of 0. We refer 
the reader to the two lecture notes of Yor [E] ) and of Mansuy and Yor |SJ , where a chapter 
is devoted to the subject. 

Inspired by some elementary remarks on renewal processes, we present an invariance 
principle for Azema martingales, and propose a new device to construct solutions of Emery's 
structure equations. The principle of these approximations turns out to be quite unusual 
and differs from the truncation method one can find in dealing with Levy processes, or from 
the discretization scheme already used in the subject by Meyer in ^U] to construct solutions 
of Emery's structure equations. It is mainly based on the introduction of some randomness 
in the size of the jumps together with keeping the normal martingale property. 

2. The basic example: the first Azema martingale 

The aim of this section is to introduce, in a natural way, the device we use in the next 
section to approximate general Azema martingales. 

Let us start with the first Azema martingale. It was obtained by Azema in by pro- 
jecting a Brownian motion Bf starting from 0, on the filtration of sign (B). This projection 
yields the following normal martingale 



Xt = sj^s^an {Bt)y/ {t - gt) 



where gt = sup{s : s < t, Bg = 0}. We refer the reader to Chapter IV of Protter's book Jl] 
for a comprehensive discussion about the martingale property of Xt- A clear justification 
of the projection property is given in an article of Azema and Yor 0: they use an elegant 
path decomposition of the Brownian motion on [0, t] involving the Brownian meander, seen 
as the renormalisation of the last incomplete excursion. 
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Now, it is well known, from P. Levy, that 2X1 jixt follows the arcsine law. On the other 
hand, there is an old result of Dynkin presented also in Feller's book Chap. XIV.3), 
obtaining arcsine laws from asymptotic waiting times of renewal processes having non- 
integrable inter- arrival times. Now, the question is: can we associate a martingale to any 
given renewal process, exactly like the Azema martingale is associated to the zeroes of the 
Brownian motion? 

The following elementary proposition gives a positive answer to this question: 

Proposition 1. Let Xi,X2, ... he an iid sequence of positive random variables admitting a 
density function. We denote by Sn ■= Yli=i -^i> '^''^^ Nt = sup{n : Sn <t}. Now, we denote 
by F{t) = P{Xi > t), the tail distribution function of the Xi 's. 

Let us now introduce a sequence ei,e2,"- of independent symmetric Bernoulli variables 
i.e. satisfying P{ei = 1) = ^ (ind P{ei = —1) = |. 

The two following assertions hold: 

(i) For every law of Xi, the process Zt ■= 7 — ^^^7; — r is a martingale with respect to its 

0F{t - SatJ 

natural filtration. 

(a) The martingale Zt is normal if and only if Xi follows the distribution of 

-lnC/+ At - - 
2[/2 2 

where U denotes a uniform variable on [0, 1]. 

Proof: The proof of (i) is based on the fact that, conditionally on the event (i — Sn^ = x), 

0F'(x) 

the probability for Zt to jump between t and t -\- h is equal to ^. . h + o(h), and if a 

0F[x) 

jump occurs, its expectation is equal to — — On the other hand, when there is no jump, 

0F 



the increment of the process, between t and t + h, is equal to —ENt t^, wo ^ + oih). As a 



0F'{x) 
{0F{x)f 

result, the expectations of the increments in both situations are balanced. 

Let us now prove (ii). The martingale is normal if and only if the conditional expectation 
of the square of the increment of Zt between t and t + h\s always equal to h + o{h). Again, 
conditionally on the event (t — Sn^ = x), the probability for Zt to jump between t and t + h 
0F'(x) 

is equal to ^, , h + o(/i), and if a jump occurs, the expectation of its square is equal to 

0F[x) 

— ((1 H )^ + (1 )^) which is 1 H ^7. Now, the absence of a jump contributes only 

2^^ 0F' ^ 0F' ' 0F2 ' ^ ^ 

in a o{h) in the expectation of the square of the increment of the process. 
Therefore, the martingale is normal if and only if 

This yields — \n0F{x) H „ = x, so that the distribution function of Xi is the 

20F^(x) 2 

inverse of x 1— > — ln(l — x) -\ — ; -75- . This gives the result. □ 

2(1 — x)^ 2 

Remark: We might have included in Proposition 1, analogous results concerning the 

vr 



sometimes called "second Azema martingale" i.e. \l —^/t — gt — k, where k denotes the 
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local time of zero at time t. Indeed, for every law of Xi, Zt := —pr-, — r- — A't is a 

0F{t - SnJ 

martingale which is normal if and only if X-i follows the distribution of ^ . 

& J 1 2C/2 2 

But we shall not insist on that fact since this second Azema martingale is neither mar- 
kovian nor enjoys the chaotic representation property, and does not enter the class of Azema 
martingales solving structure equations. 

Let us denote by zj: the normal martingale Zt of (ii). We note that zj: is a Markov 
process, and more precisely. 

Proposition 2. The process z\ is a Markov process with generator: 

Ug{-l)+gil))-g{x)+xg'{x) 



L(-i)<7(x) 



1 + x2 



Proof: From Proposition 1 (ii), we deduce that 



-ln0F{x) H ^„ , , = X 

^ ' 20F2(x) 2 



which implies, by differentiation: 

0F'{x) 



0F{x) I 

0F2(x) 

Hence, the rate of jump of the process at time t, which is equal to the value of the function 
0F' 

^■t (t ~ ^Nt): is precisely ^ ^(-i)-^2 ' ^^'^ speed of the trajectory at a time t 

0F'(t-SNt) . , , ^^^^ 

between two jumps is equal to —SNt — 7, — r which is precisely i—r: — . □ 

0F^{t-SNt) l + {Zl 'f 

(-In) z'-~^^ 

Corollary 1. The process Z\ := ^ is a Markov process with generator: 

U9{-j^) + 9{j^))-9(x)+xg'{x) 



Proof: Consider a differentiable function g. Introduce the function h{x) = g[-^). 

The image by L^"^'") of the function g is given by: 

L(-i'")5 = n.(L(-i)/i)(V^.x) 

A direct computation gives the announced expression. □ 

This last corollary already gives an idea of the proximity between Z" and the first Azema 
martingale. Indeed, as writes Emery in |5j about Azema martingales: "formally (this means: 
informally !) it should be Markov, with generator 

r g{{l + l3)x)-g{x)-l3xg'{x) 



Lg{x) 



22 if X / 

y{x) ifx = o" 



I3^x^ 

1 J' ( 



(in the case of the first Azema martingale, the parameter (3 is equal to —1) 

This is the aim of the next section which is to prove, among other things, the convergence 

^(-1) 

in distribution in the sense of the Skorohod topology of ^ towards the first Azema 
martingale. 
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3. An invariance principle 

Let us begin with the definition of the Azema martingales. 

Following Emery, let us introduce a real parameter /3. The structure equation 
(1) d[X,X]t = dt + pXt-dXt 

has a unique solution. 

The existence was proven in a note of Meyer jlUj . and the uniqueness of the solution 
was proven by Emery in :5_ (for /? < 0) and [H] (for /? > 0). The solution of this structure 
equation is called the Azema martingale with parameter /?. 

When a jump occurs at time t, the value of the process changes by some prescribed factor. 
More precisely, Xt = {1 + (3)Xi-^ . 

• The case /? = corresponds to Brownian motion: the continuous normal martingale. 



• The case /? = — 1 is solved by the martingale {sign Bt)\/2{t — gt) mentioned in the 
previous section. 

• The case f5 = —2 is solved by Parthasarathy's martingale satisfying \Xt\ = \ft^ and who 
changes of sign according to a Poisson point process with intensity dt/4t. 

• In the case (3 < —2, the solution, at least starting at a non-null point, is rather easy to 
define, since \Xt\ always goes away from 0. In proving the uniqueness of the solution, Emery 
in j^l shows how to define the process starting from by using a self-similarity argument 
and a representation of the process by using a time-changed Poisson process. 

For P > —2, things go in a more complicated way since the process Xt reaches in finite 
time, and the above formal generator cannot be of much help after that time. 

Let us mention finally that the problem of the chaotic representation property still remains 
a challenging open question for parameters (5 which do not belong to [—2, 0]. 

Now, let us turn to our main purpose concerning the statement of an invariance principle 
for these processes. This question is quite natural, since Azema martingales are self-similar 
in the sense that, for all A > the processes Xt and Xx2t/\ are equal in distribution (one 
can indeed easily check that the process X^it/X solves the right structure equation, and 
conclude by uniqueness in law of the solution). 

In order to define the process which will be, in our statement, at the origin of such a limit 
theorem, we generalize the definition of zj: introduced in the previous section, and define 
for all real parameters f3 the Markov process zj:^^ having generator: 

^ ^(g((l + f3)x - 1) + g((l + (3)x + 1)) - g{x) - (5xg'{x) 

1 + /?2x2 

(B n) Z*"' 

Let us denote the process Zi := It has generator: 

_ 1 (g((l + (3)x - ^) + g((l + (3)x + ^)) - g{x) - Pxg'jx) 

^ ' " ^ + f3^x^ 

We can now state our main result: 

(0 n) Z'-'^^ 

Theorem 1. The sequence of processes Z^ ' := converges to the Azema martingale 
with parameter (3, in the sense of the weak convergence for the Skorohod topology. 

Proof: The image by L^^'^\ oi x and x^, being respectively and 1, we deduce that zj.^'""^ 
are normal martingales. Therefore, by Rebolledo's theorem (see jl5| . II.3.1), this sequence 
is tight for the weak convergence in the Skorohod topology. All we have to prove is that all 
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limiting points of the sequence of processes Z^^'^^ satisfy the structure equation (1). The 
uniqueness in law of the solution of (1) allows to conclude. 

Let us define at a jump time t of the process Zf^'^\ the symmetric Bernoulli variable 
s[^'"''' defined by 



(/3,n) 



'n 

Let us consider [Z^^'^\ Z^^'"')]^. (We replace, in the sequel, the superscript (/?,n) by n.) 
[Z-, Z% = Y,{Z^ - Z:-f = + -j=){Z- - 

s<t s<t 

= / pZ-.dZ:- (3Z-.l^z^^odZ2 + Y.^Azn^<^^{pz:. + ^) 
Jo Jo 

t r-t J pTb 

pz';.dz: +t- ' + lAz^^o^iPz:. + ^) 

Jo l + n{f3Z^_y ^ ^/n ^Jn 

Our task, now, is to show that the two last terms vanish when n goes to infinity, whereas 
the other terms converge to their analogous quantity for the limiting process, so that every 
limiting process will solve the structure equation satisfied by an Azema martingale, and will 
therefore be identified with the Azema martingale with parameter /3. 

Now we proceed, like in Meyer's proof of the existence of solutions of structure equations, 
and introduce successively: 

_ a probability space (0,^,(5), obtained from the Skorohod representation theorem, on 
which a subsequence of processes Xl^"'\ following the law of an extracted sequence of 
converges almost surely in \ ®Q. For ease of reading, we keep writing instead of Xl^^^ . 
In other words, we suppose that, for almost every time Z^ converges almost surely. 

_ a sequence of stopping times = inf{t : > A''} where N is large enough, so that 
Q(T^ < is smaller than any prescribed e, uniformly in n (Doob's inequality makes it 
possible to do so, since the processes have the same variance). 

Let us compare, in a first step, [Z", Z^jxj^^t with the sum of the square of its increments 
along some subdivision (tj) of the interval [0,t]: 

i i 

The compensator of Ztj^as being equal to A s, we deduce, 

EmZ-, Z-]TnAt - E(^?;^At,+l - Z^.^,ff] < 4^[E (^Tn A._ - Z^n^,fds] 

which in turn is smaller that 'iElY,^ //^'+' (^r;jAs_ " ^r^jAtJ^^^] ^ I T^iiU+i - Uf- (In this 
first step, stopping the martingales at just ensures the integrability we need. In the 
further steps, this stopping argument will be more widely used.) 

As a conclusion, we deduce that the "discrete quadratic Riemann sums" of Z" approx- 
imate [Z", Z"] uniformly, in probability. Now, since the corresponding sums for Z ap- 
proximate [Z, Z] in probability, it just remains to notice that we can choose the time of 
our subdivisions among the times t for which Z^ converges almost surely to Zj so that 
the quadratic sums of Z" will converge almost surely and therefore in probability to their 
corresponding analog for the process Z. 
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In a second step, we deal with the second term Jq Z^_dZ^. We prove its convergence to 
the corresponding quantity for Z, exactly the same way as in Meyer's proof. 

Now, we have to prove the convergence to in probability of both remaining terms 

ft j„ pU _n 

Jq l + n{(3Zl'_y ^ Vn ^/n 

Again, it suffices to prove it for these quantities when t is replaced by A t. 

This result will be a consequence of the following key lemma, which gives some uniform 
control (in n) on the time spent by the processes near the origin, and on the number of 
their jumps: 

Lemma 1. The two following convergence results hold: for all N > 0, and all t > 0, 
(i) Ve > 0, 3(5 > 0, 3no G N, such that Vn > no, 

P{ I lzr.^\_f,ji\ds > e) < e 

^0 



LZ,"e[-<5,<5] 



T" At 

(a) — - — converges to in probability, 
n 

Proof: We start with (i). We will use the following notation for a constant which will 
appear repeatedly: C = max(l, 1 + 

Let us fix e > 0, and, for some 6 g]0, ^[, consider the successive stopping times: 

_rf = inf{t : > 25} 

_ = mi{t > rf : |Zf | < 6 or \Zf \ > N} 

_ If T^' ^TJ^, wc define t^' := mf{t > : \Z't\ > 26} 

and more generally for i > 2, 

.T^,:=mi{t>T^,_,:\Zl^\<6 or \Z,-\ > N} 

_ If T^^ + T^, then r2Vi := mf{t > r^', : i^fl > 2b} 

Let us denote by the random integer K := inf{i > 1 : r^^j = T^}- 

Prom these definitions, we get obviously: 

r^'^' lzne[-5,5]ds < rr + {r^ - r^) + ... + {t^k-i - 
Jo 

Now, we can bound the expectation of the right hand term, by first noticing that, by the 
optional sampling theorem applied to the martingale {Z't)'^ — t, 

Vi>l, E[r^,^,-TS,\K>i] =E[{Zl]^^J^-{Zj;^f\K>i] 

< E[{Z';n^^f \K>i]< 4max(l, (1 + pf)6^ = 4C^d^ 

And secondly, we can estimate P{K = i + 1 \ K > i), hy exploiting that 

ElZj^n \K>i]= E\Z\ \K>i]. 

'2i+2 ' ' '2i+l ' ' 

But, on one hand, \Z'!ln I > 25 and, on the other hand, ElZJln \K > i] decomposes 
into E[lK=i+iZ'^n ]+E[lK>i+iZ^n ],with 

\E[lK=i+iZ%]\ < NC.P{K = i + l\K>i) 



'2i+2 

and 



\E[lK>i+iZ:^^^J\<5 
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Therefore, 

r 

P{K = i + l\K >i)> 



NC 

which impUes that P{K > i) < {1 - j^Y^^- 
Now, writing 

i>i 

and, using the strong Markov property at time T^j_i, we get 



i>i 

_ ^2 



We just have to choose S such that 4NC S = e and use Markov inequahty, to get (i). 

To prove (ii), we will strongly exploit (i). To make the reading easier, we will fix t equal 
to 1. Let us divide the interval [0, t] = [0, 1] into n intervals of the type -^^ = ^^], where 
< A; < n - 1. 

Now, e being fixed, we introduce 6 like in (i), and define the stopping times: 

= inf{. G 4", \Z^\>S} 

with CT^ = 00 if Vs G \Z^\ < 6. 

From (i), we get that, for n large enough, with probability bigger than 1 — e, for more 
than (1 — e)n integers A;, between and n — 1, we have crj^ < 00. 

Let us introduce now the stopping time := inf{s > cr^ ~ ^s- 1 7^ 0}) "^^i^ '^k — ^ 
whenever cr^ = 00. 

Now, for n large enough, the quantity P{\tJ^ — < - | o'fc < 00) can be bounded 
uniformly in k by some sequence depending on n, which converges to 0. Indeed, if at some 
time s, \Zg \ > 5, and if moreover there is no jump between s and s + f , then \Z^\ remains 
bigger than 5/2 for n large enough. Now, on the time interval [s, s + ^], the intensity of the 
jump remains smaller than ^t^jj^- 

We deduce that, conditionally on the event (a-^ < 00), — (^k \ ^ n stochastically 
bigger than 6 4^ A ^, where 6 4^ is an exponential variable having parameter pTg^- 

So that, for n large enough, 

2 8 
^(kfe - (T^l > - I < 00) > exp(--^2^) 

Now, — (T^l > ^ implies that there is no jump inside the time interval 

So, we deduce that, for n large enough, with probability bigger than 1 — 2e, the process 
Z^ will not jump on more than (1 — 2e)n intervals among the /^'s. 

Finally, since the rate jump of the process is permanently less than n and the length of 
each /" is equal to ^, for n large enough, we get that, with an arbitrarily high probability, 
there are not more than 4en jumps on the remaining 2en intervals which may have some 
jumps. □ 

End of the proof of Theorem 1: The (i) part of Lemma 1 gives immediately the convergence 

m probability to 0, of / — — , 

7o l + n(/3Z!»_)^ 
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Now, V u^„^o>(/3^r- + >) = ^^+/3 E lAz^vo^^r- 

.<^M ^ ' ^ n ^^^^^ ^ 

Part (i) of Lemma 1 implies the convergence to of the first term of this sum, and we 
bound the L^-norm of the second term: 

s<TJ^At s<TJ^At 

Now, for all e > 0, 

ATn ivrn njn nrn 

E[-^^] < e + Ei^i^.l ^n^^ ]<e + E[{^^f]-2 Pi^^ > e)^ 



Now, coupling the jump times of the process with the points of a Poisson process with 

AT" 

intensity n, we are able to bound E[{ — ^ — ) ] uniformly in n, by t . 



We conclude that 1az"^o—7={PZ"'- + —^) converges to in probability, and obtain 

Theorem 1. □ 
Remark: We might have replaced the randomness of the jumps, having distribution 
2(^-1 + 61) in the case of the process , by any other one, having compact supported 

distribution with mean and variance 1. In other words, the result of Theorem 1 remains 
valid, if we take, instead of Zf^'^\ any process having generator 

+ P)x + y)fi{dy) ] - g{x) - f3xg'{x) 



1 + /32,x2 

where is a compact supported measure, with mean and variance 1. Indeed, one can 
easily check that we still deal with normal martingales. 



4. An extension to structure equations 



Actually, Meyer is able to construct solutions of the more general structure equation 

(2) d[X,X]t = dt + f{Xt-)dXt 

where / is an arbitrary continuous function. 

The strategy of Meyer consists in solving a discrctizcd structure equation, which leads 
to a cascade of second degree polynomial equations satisfied by the discretized increments, 
having always exactly two solutions of opposite signs, which are chosen with respect to the 
(unique) probability preserving the martingale property of the constructed discrete process. 

We want to show in this section how our approach can be used to exhibit a solution of 
the structure equation (which may be different from Meyer's in the non- uniqueness cases). 
Inspired by the processes of previous section, we introduce a sequence of normal martingales 
zl^'^^ whose limiting points in the sense of the weak convergence for the Skorohod topology 
will satisfy the structure equation (2). 

However, although our approach has the advantage to give an invariancc principle in 
the self-similar case of previous section, nonetheless the construction of structure equations 
for general continuous functions induces delicate discussions on the time spent near the 
zeroes of the function /, and wc did not work them out for a general continuous function 
/. So, our aim, in this section, is more to examplify things, than to arrive to the most 
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general statement, which might be vahd for general continuous functions, but we will let 
this question open here. 

We define zj:'^'"'\ as the Markov process starting at 0, and having generator: 



Proposition 3. Suppose f is a continuous function with isolated zeroes. Suppose, in addi- 
tion, that, in the neighbourhood of each zero Xj of f, f{xj + h) = o(\//i)- 

Then, the limiting points of the sequence of processes zj:^'"'\ in the sense of the weak 
convergence for the Skorohod topology, satisfy the structure equation (2). 

Proof: Like in the previous section, we start by noticing that the image by L^^'"'\ of x 
and x^, is respectively and 1. Indeed, 

1 ((x + fix) + ^) + (x + fix) - -^)) -X- fix) = 

2 V \/n Jn / 



and 

\ [ix + fix) + + (x + fix) - - x2 - 2x/(x) = i + fixf 

2 V Jn Jn / n 



We deduce that zj:'^'"'^ are normal martingales. Therefore, by Rebolledo's theorem (see |15j . 
II. 3.1), this sequence is tight for the weak convergence in the Skorohod topology. 

Now, we follow the proof of Theorem 1, and write: 

ft r-t J pTi 

Applying again Skorohod representation theorem, we are able to prove the convergence of 
the terms [Z", and fiZ^_)dZg to the analogous quantities for the limiting process. 

Now, the key point is to check the analogous property to Part (i) of Lemma 1. 

In other words, we have to prove, that, for every xj G [— A^, N] (there are a finite number 
of them), Ve > 0, 36 j > 0, 3nQ £ N, such that Vn > no. 



P{ / lz^^e[x,-5,,x,+5,]ds > e) <e 
Jo 



We then introduce the same stopping times as in Lemma 1, and the differences with the 
proof of Lemma 1 come out when we estimate -EiT^j+i — T"2i \ K > i]: 

yi > 1, E[t^^^, - T^, \K>i] = E[iZ!^^^J^ - iZ^^^f \K>i\ 

< EiiZ^^^J^ \K>i]< (25, + sup[„25„25,] 1/(^)1)' = oi5j) 

whereas PiK = i + 1\K > i) remains bigger than a constant times 6j. 

As a result, E[ti + (r^* — r^^) + ... + (t^/^_-^ — t^;^_2)] is bounded by a constant times 

(2<5j+sup[_2i 2^,-1 1/(2^)1)^ 
^ = o(l)- 

As in Lemma 1, we choose 6i such that this fraction is equal to e^. 

We finish the proof the same way as for Theorem 1. □ 

Remark 1: The framework of Proposition 3 includes the case of asymmetric Azema mar- 
tingales corresponding to a function / of the form axla;>o + bxlx<o- In this case, Phan 
proved in [12j that the solution of the structure equation is unique. From this result, we 
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deduce that these processes are self-similar and Proposition 3 provides, in this case, an 
invar iance principle. 

Remark 2: Our assumptions in Proposition 3 are obviously not optimal: for instance, 
if / is null, we are back to the case of a simple random walk, and the convergence to the 
corresponding solution of the structure equation which is, in this case, the Brownian motion 
is obviously valid. It would not be hard to check that one can allow / to be null on a locally 
finite number of intervals, and satisfy the assumptions of Proposition 3 on the complement 
of the union of these intervals. In this case, the term X]s<t lAZy7^0"^(/(-^"-) + will 

i /(Z" )^ 

not vanish in the limit and will "alternatively" contribute with j_ ^~ ^ ds to build a 

compensator equal to t. 

It is also very likely that the condition on / to be o(v^) near its zeroes is not optimal 
(even if our attempts to relax it failed). 

5. Construction of the process z^^^'") 

We might just present the construction of zj:^'^-' and rescale it, but we prefered to keep 
the parameter n in evidence, in order to point out clearly where the randomization of the 
jumps appears. 

Let us introduce the functions 0F and H defined on M x R+, in the following way: for 

any (xo,t) G M x R+, 

_ 0F(xo,t) denotes the probability that the first jump of the process Z^^'"'\ starting at 
xq, takes place after time t. 

_ H{xo,t) is the value of zj-^'"'\ starting at xq, conditioned on the event that the first 
jump of the process takes place after time t. 

Prom the expression of the generator, we deduce the following system of differential 
equations: 

1 d0F 1 
-{xo,t) = - — 



0F df"' ^ + ^m^ixo,t) 



with the boundary conditions 0F(xo, 0) = 1 and H{xq, 0) = xq. 

. . , 1 d\nH d\n0F , , 
After noticing that — — - — = — , we deduce that 

p ot ot 

H{xo,t) = xo{0Fixo,t)f 
Prom this fact, we deduce the following autonomous equation for 0F: 

1 1 d0F ^2 2 ^2fl ^d0F ^ 
^-7^ P X^0F^f^-^-^ = 1 

n0F dt ^ ° dt 

which yields 

_lln0F-^(0F2/3-l)=i 
n 2 

Therefore, the first jump of the process zj:^'"'\ starting at xq, follows the law of 

_iln[/-M([/2/3_i) 
n 2 ^ ^ 
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where [/ is a uniform random variable on [0, 1]. We notice here that the special case n = 1, 
P = —1 and \xo\ = 1 allows us to recover the law of the random variable we discussed in 
Proposition 1 (ii). 

Denote Tq = and Tj the i-th time of jump of the process We can describe the 

process zf'^^ as follows: there exists a sequence (J7j)j>i of independent uniform variables 
on [0,1], and a sequence {£i)i>i of independent symmetric Bernoulli variables, such that, 
for all i > 1, (we omit here the superscript (/?, n).) 

( 1 z'^ 

Ti = Ti_i--lnU,-P^iUf^ -1) 
< Z^^ = ZT,_,uf 

1 

Zt^ = Zy-(l + /3) + — 

On the time interval [Tj_i, Tj[, the trajectory (t, Zt) can be parametrized by -^((s))o<s<i-(7i 
in the following way: for all s S [0, 1 — Ui[, 

1 z"^ 

t{s) = r,_i--ln(l-.)-/3^((l-s)2/3-l) 

Zti^s) = ZT,^,{l-sf 

Remark: From this description of z[^'^\ we can easily deduce a coupling between z[ ^'"^ 

and z[^'^\ for all n and all parameters /?. But, unfortunately, this coupling cannot be 
transported at the level of the corresponding Azema martingales by merely making n tend 
to infinity. 

This description allows us to draw some pictures of the process zj:^'^^ for large values of 
n and different values of the parameter (3, having therefore some good approximations of 
trajectories of Azema martingales. However, in this paper, we do not intend to quantify 
any error term. 

We finish this paper with some simulations of Azema martingales, which were obtained 
by writing a code in MATLAB. It basically draws z[^'^\ asks for the number of arches one 
wants to appear in the graph, and rescale the obtained function in order to get a function 
defined on [0, 1]. In the figures below, we asked for 1000 arches to appear. 

Figure 1 represents Parthasarathy's martingale, whose graph is included inside the parabola 
of equation \y\ = ^fx. Figure 2 represents the first Azema martingale. We can see on Figure 
4 that, for very small values of the parameter /3, the Azema martingale is a perturbation 
of the Brownian motion. It was considered as a candidate in ^ to replace the Brownian 
motion in the Black-Scholes model. We chose /3 = —0.07 which was the parameter formerly 
estimated by Dritschel and Protter to modelise the variations of some asset prices. 

We finally couple figures 3 and 5, which correspond respectively to the parameters —0.5 
and 1, and notice similarities in their graphs up to some symmetry and horizontal affinity. 
Parthasarathy conjectured, after noticing the duality between the two formal generators of 
these processes that, up to a time change, the first process was the time reversal of the other 
(the conjecture was formulated, of course, for general positive values of the parameter). This 
conjecture was actually proven in an unpublished part of Phan's thesis jj^ by a nice and 
unusual argument involving polynomial test functions. 



N. ENRIQUEZ 




Fig. 4: /3 = -0.07 
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Fig. 5: (3 = 1 
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